Unstable fingering patterns of Hele-Shaw flows as a dispersionless limit of the KdV 
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We show that unstable fingering patterns of two dimensional flows of viscous fluids with open 
boundary are described by a dispersionless limit of the KdV hierarchy. In this framework, the 
fingering instability is linked to a known instability leading to regularized shock solutions for 
nonlinear waves, in dispersive media. The integrable structure of the flow suggests a dispersive 
regularization of the finite-time singularities. 
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1. Introduction. The Hele-Shaw cell is a narrow 
space between two parallel planes - a device used to study 
the 2D dynamics of a fluid with an open boundary. In a 
common set-up, air (regarded as a non- viscous fluid) oc- 
cupies a bounded domain of the cell, otherwise filled by 
a viscous incompressible fluid. When more air is injected 
through a well, the free boundary evolves in a compli- 
cated and unstable manner. In finite time, an arbitrary 
smooth initial boundary develops a pattern of branched 
fingers . This mechanism is rather general and has been 
identified in numerous growth problems Q . 

Starting from early works 0] , it became clear that the 
problem of fingering instability in the Hele-Shaw cell is 
linked to profound aspects of the analytical functions the- 
ory and 2D conformal maps. In recent years, the prob- 
lem has reached other domains of theoretical physics. In 
particular, a connection to fractal self-similar patterns of 
stochastic growth, like DLA, has been recognized Q- 

The Hele-Shaw problem has also emerged in the con- 
text of electronic physics in low dimensions. For ID elec- 
trons, the Fermi sphere is deformed by externalpertur- 
bations according to the Hele-Shaw mechanism |4j. The 
same is true for deformations of an electronic droplet in 
a quantizing magnetic field ||. Similarly, the Hele-Shaw 
problem is connected to statistical ensembles of normal 
or complex random matrices, where a non-Gaussian po- 
tential deforms the support of eigenvalues jg]. 

In this paper we emphasize yet another connection, 
which may be viewed as a formal ground for appearance 
of this phenomenon in different physical situations. We 
show that the fingering instability is linked to the inte- 
grable Korteweg-de Vries (KdV) hierarchy of differential 
equations. More precisely, we show that developed fin- 
gers, in the absence of surface tension, are described by 
the dispersionless KdV hierarchy. The Hele-Shaw fingers 
correspond to the same solution of the KdV hierarchy as 
that exploited in 2D-gravity , and in studies of critical 
points of Hermitian random matrix ensembles 8] . 

The integrable structure of the Hele-Shaw problem has 
been observed in our earlier papers Q. There, we have 
shown that the problem is equivalent to the dispersionless 



limit of the Toda integrable hierarchy. This is an exact 
relation. Here we concentrate on the critical (turbulent) 
regime of the flow, where fingers are already well devel- 
oped, and close to a cusp-like singularity. In this case, 
one may concentrate only on the vicinity of the tip of a 
single finger, neglecting the rest of the boundary. The 
KdV integrable structure emerges in this regime. It can 
be obtained from the Toda hierarchy as a reduction, but 
can also be derived directly from hydrodynamics of the 
critical regime. We choose the latter. 

In this letter, the intuitive physics of interface dynam- 
ics gives a clear geometrical interpretation to algebraic 
objects used in dispersionless soliton equations [lfj. 

2. Darcy's law. Consider a 2D domain (a "bubble") 
occupied by "air" regarded as an incompressible fluid 
with low viscosity. The rest of the cell is occupied by 
another incompressible fluid but with high viscosity. Air 
is injected into the cell through the origin, at a constant 
rate, while the viscous fluid is evacuated from "infinity" 
(edges of the cell). The area of the bubble is proportional 
to time t. We normalize it to be nt. The Navier-Stokes 
equation adapted to the 2D geometry gives a simple rule 
for the dynamics of the moving boundary 7: velocity in 
the viscous fluid (and of the boundary) is proportional 
to the gradient of pressure, 



v = — Vp. 



(1) 



Pressure is harmonic for incompressible fluids. Inside the 
bubble, it is constant (set to p = 0), due to low viscosity. 
In the absence of surface tension, it is continuous across 
the boundary, and hence a solution of the problem 



Ap = 0, p\ 



log I 



P\ = 0. 



(2) 



The origin of the fingering instability is intuitively clear 
- high curvature portions of the boundary move with a 
higher velocity than the rest and get even more curved. 

We now introduce a minimal set of objects of poten- 
tial theory. Consider a holomorphic function 4>{z) = 
£(2, y) + ip(x, y), whose imaginary part is pressure. The 
function <f>(z) is a univalent conformal map of the exterior 
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of the bubble to the cylinder Im<p > 0, and its derivative 
d<p(z) taken on the boundary is the conformal measure of 
the boundary. Darcy's law Q reads: the complex veloc- 
ity of the boundary is proportional to the conformal mea- 
sure. In other words, the complex potential of the flow 
is a conformal map of the outer domain to the cylinder. 
Equivalcntly, the Cauchy-Riemann conjugate of pressure, 
£(cc,y), is the stream function. 

For simplicity, let the bubble and finger have a sym- 
metry axis. In Cartesian coordinates, the boundary can 
be described by a multivalued function y(x), Fig^ 
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FIG. 1: A finger extends to the right of its boundary 
y(x). —2u(t) is the distance from the tip to the singular point. 



Another useful way to describe the boundary is 
through the height function h(z). It is an analytic func- 
tion in some vicinity of the boundary, taking values 
h(x + iy) = y(x) on the boundary . With the help of the 
Schwarz reflection principle, we obtain h(4>) — h(—(f>) = 
—i(z((f>) — z(— </>)), where z((f>) is the inverse of the map 
(we write h((f>) — h(z((f>))). Therefore, 



2ih{4>) = *(<t>)-z(- 



(3) 



We will also use a generating function d z uj{z) — ~2ih(z). 
Darcy's law can be interpreted as evolution of the height 
function. The complex velocity v = z can be analyti- 
cally continued into the outer domain, as v — dtd z uj(z) — 
—2idth(z). Then from v = —2id z <f), Darcy's law becomes 



d t h(z) = d z (j)(z). 



(4) 



3. Branch points If the boundary is analytic, the map 
4>(z) and the height function h(z) can be analytically 
continued inside the bubble, until they reach singulari- 
ties. For an important class of analytic boundaries, the 
conformal map z(4>) is a rational function of . Such do- 
mains are called algebraic. In this case, the singularities 
of the inverse map are generically simple branch points 
2o (£),... , z g _i(i) (analogs of Riemann invariants in hy- 
drodynamics), where d z (j) ~ (z — Zi(i))^ 1 / 2 . Hele-Shaw 
flows preserve the class of algebraic domains 11]. 

As the boundary moves, so do the branch points. We 
identify the traces of the branch points with the branch 



cuts. Comparing singularities of we get a system of 
coupled equations for the branch points, 



dh 



i = 0,l,...,q-l. 



(5) 



4- Finite time singularities In full generality, these 
equations are difficult to analyze. Qualitatively, branch 
points tend to move towards the boundary. As a result, 
the boundary is pushed away, forming a finger. When, 
eventually, a branch point reaches the boundary, the fin- 
ger forms a cusp-like singularity. It appears that for al- 
most all initial algebraic domains, a cusp-like singularity 
emerges at a finite time, i.e. at a finite area of the bubble. 
The idealized law Q then needs to be corrected by reg- 
ularizing the singularity ,12|. Some examples of singular 
Hele-Shaw flows were elaborated in |l3j. 

In this letter we suggest a unified approach to the fi- 
nite time singularities based on a singular limit of the 
KdV integrable hierarchy. We concentrate on the criti- 
cal regime of the flow, when fingers are close to a cusp. 
Generically, an isolated branch point, say zq — (xq, 0), is 
found very close to the boundary (tip of the finger), i.e., 
4>{ z o) ~ 0- From (JjjJ, this branch point (and tip) moves 
with velocity xq = d<j>/dy\ XQ . 

The origin of the finite time singularities can already 
be seen from this equation by means of scaling analysis. 
Around the tip, the map z(<j)) has a regular expansion 



z(cf>) — xq = iacj) + b(j)- + icifi + 



with real a,b,c . 



It follows from the reflection formula © that h(cf>) = 
a<f> + ccj) 3 + . . . . If a branch point is close to the tip, then 
a is small, while b and c are of order 1. If there is only one 
scale, u(t), then xq, measured from what will be the cusp 
tip, and a are of the same order u(t). The velocity of the 
finger tip xq ~ u grows with the curvature a~ l ~ u^ 1 of 
the tip such that u ~ u~ l , hence the scale vanishes as 
u ~ (t c — t) 1 / 2 . At the critical time t — t c , u = and the 
curve forms a (2,3)-cusp, y(x) ~ x 3 / 2 . 

More general cusps are characterized by two integers 
(q,p), q < p, implying that the finger is bounded by 
\y(x)\ ~ x p / q , x > 0, q is even. Higher cusps correspond 
to higher order singularities merging at the boundary. 

In what follows we mostly consider the case q = 2, 
making only brief remarks on the higher cusps. The de- 
tails will be given elsewhere. 

5. Critical regime and scaling functions. When the 
finger is very close to becoming a cusp, the rest of the 
bubble does not affect its evolution. We will call this 
regime critical. In the critical regime, different scales 
separate. It is convenient to define the scales in the com- 
plex plane of <f>. At cf> = (^(u 1 / 2 ) we see the details of the 
tip. At 1 3> 4> ^S> u 1 / 2 , the details of the tip cannot be 
seen, the finger looks like a cusp. Finally, <f> = 0(1) cor- 
responds to the rest of the body of the bubble. The scale 
u(t) changes with time and eventually disappears at the 
critical point. In that limit, the finger is scale invariant. 
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To summarize, in the critical regime u — > 0, t — > t c , 
the inverse map and the height function are 

z{(t))-ih{(j}) = u^Q^/u 1 ' 2 ), h{(j>) = u^P{<j)/u 1 ' 2 ) (6) 

where the scaling functions Q and P do not depend on 
u. At q — 2, Q is an even quadratic polynomial and the 
scaling ansatz holds everywhere in the critical region. 

The asymptotic behavior of the finger, y ~ ±x p / q , al- 
lows us to identify the scaling functions. They are poly- 
nomials of degrees q and p, respectively. The reflection 
symmetry suggests that Q is an even polynomial, while 
P is odd, and that all their coefficients are real. This 
follows by matching asymptotes in different regions: far 
from the tip, the finger asymptote is y ~ ±x p / q . Since 
j/CiCl, we approximate y ~ z p / q , where z = x + iy. 
Higher order corrections will make y real. Therefore, we 
conclude that the height function, in the domain close 
to the boundary, but away from the tip, is h(z) ~ z p ^ q . 
This also means that at u 1 / 2 <C 4> "C 1, and close to 
the real axis, Q{4>) ~ </> 9 , P{4>) ~ </> p - Close to the tip, 
where <C w 1 / 2 , both z(^>) and /i(<^>) are regular in <f). 
The only holomorphic functions matching these condi- 
tions are polynomials. 

6. An important remark on the scaling ansatz is 
in order. From 0, approximating the univalent map 
z (</)), being by polynomials makes it no longer univalent 
at q > 2. It covers the z-plane q times, as many as 
the number of sheets of the complex curve. This prob- 
ably means that at q > 2 the scaling ansatz does not 
hold everywhere in the upper half-plane of ef>. It holds 
only in certain sectors of the plane adjacent to the real 
axis and breaks down otherwise. Their union, covering 
roughly 2/q part of the upper half-plane - a physical 
domain determines a physical branch of the map 4>{z). 
The physical branch must obey the reflection symmetry 
and be real at the boundary of the finger. This is impos- 
sible unless the stream function Re<j) is allowed to have 
a discontinuity across the reflection axis outside the fin- 
ger. Every monomial z n l q in the asymptotic expansion 
4>{z) — z x l q + negative powers of z x l q should be under- 
stood as z 1 ^, if < Argz < tt - e, and e 2vi / q z x / q , if 
— 7r + e < Argz < 0. As a result the imaginary part of 
(j>(z) has a finite discontinuity on the reflection line, at 
least far away from the tip, unless q = 2. 

This indicates that cusps with q > 2 do not appear for 
algebraic simply-connected domains. We do not have an 
interpretation of this phenomenon. The non-univalidness 
of the map suggests that there is another finger at the 
left, while the cusp singularity implies their merging. 

7. dKdV hierarchy From 10, in the scaling limit, the 
generating function has the asymptotic expansion u>(z) — 
uj(z Q ) = -2i J* o h(z')dz' + 0{u p ) « u^ J* PdQ. It is a 
truncated Laurent series in z 1 ^, 

p+q 

u)(z) = i t n z n l q + negative powers of z 1 ^. (7) 



We will see in a moment that all the coefficients tk ex- 
cept t\ are conserved. Moreover, the coefficient t\ is pro- 
portional to time, measured from the moment of singu- 
larity. The coefficients tk-,t\ are called deformation and 
evolution parameters, respectively. The generating func- 
tion and the height function can be expressed through 
deformation parameters. As polynomials of <f), they are 

p+q n 

t n u) n , h= —t n ui n - q . (8) 



p+q 



n=l 



71=1 + 9 



Here w„(</») is the polynomial part of u n / 2 Q n / q ((j)). 

Now we can describe the evolution of the curve with re- 
spect to all the deformation parameters. The arguments 
are borrowed from 01 with little changes. We note that 
dt n w is a polynomial of degree n in <p. As a Laurent series 
in z 1 / 9 , it has only one term, z n / q , with positive degree. 
The only polynomial of this kind is co n . Therefore, 

dt n h(z) = d z ui n (z). (9) 

Notice that the flow in real time appears on the same 
footing as flows with respect to the deformation param- 
eters. Setting n = 1, we recover the flow equation J3J. 

Eqs. H4I9JI can be cast in the form of flow equations, 
if one takes time derivative at constant <p. Defining the 
Poisson bracket with respect to the canonical pair t, <p 
such that {/, g} — dtfd^g — d t gd^,f, eqs. (|4I9(I read 
{h, z} = 1, dt n h((j)) = {h, ui n }. From this we conclude 
that the inverse map evolves as 



9t n z = {z, u n }. 



(10) 



Compatibility of these equations gives a closed set of non- 
linear equations for the coefficients of the polynomials (JSJ) 
dt m uj n — dt n u>m — {wn, ^>m) ■ At q = 2, it is the disper- 
sionless KdV hierarchy (dKdV). At q > 2, equations ltTT)| 
constitute the dispersionless Gelfand-Dikii hierarchy [1J| . 

8. Solutions of the dKdV hierarchy At q — 2, a com- 
plete solution is available. In this case z(<f>) = <\> 2 — 2u, 
the evolution does not depend on tm and 



£ (2n+l)!! {-u)\ 2n _ 2k+1 



k=0 



(2n-2fc + l)!! fc! 



11) 



Equations i|10|) become 

_ (2n + l)!!. 

ot 2n+1 u = — {-u) n d tl u. 

The first equation is the familiar Hopf-Burgers equation 
d t3 u + 3ud tl u = 0. (12) 
The hodograph transformation gives a general solution 



E 

k=0 



(2fc+l)!! 

jfei 



*2fc+i(-«) fc = 0, p = 2/ + l. (13) 



More details of the application to the singular limit of 
the Hele-Shaw flow can be found in [lj|. Some results 
for q > 2 are discussed in [lij . 
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9. KdV hierarchy Eq. (|12[) is a singular h 
of the full dispersive KdV equation 



4a. 



12u<9 t u + n 5/2 u ttt = 0. 



limit 



(14) 



To clarify the nature of this limit, we recall the definition 
of the q-reduced KP hierarchy. It is a set of nonlinear 
equations compactly written through a pair of operators 
I/,M, differential operators in time of degrees q and p: 



L 



- Mi 



i 

E 

1=1 



M = 



p+q 

E : 

n— 1+q 



(15) 



Here [, ] + denotes the anti-commutator, <f> — hdt, and 
the coefficients ei are functions of a string of "times" 



. ,t p + g , and fl n — L']j_ q is a differential operator of 
degree n, obtained from the pseudo-differential operator 
L n/q by omitting a non-differential part. The dependence 
of the coefficients e; on the "times" is introduced by the 
Lax-Sato equations hd tk L = [L, Q n ], hd tk M = [M, fl n ], 
supplemented by the condition [L, M] = H. An alterna- 
tive way to define the hierarchy is to impose the consis- 
tency conditions d tn fl m — d tm fl n — [fl m , f2 n ]. 

The dispersionless ( "quasiclassical" ) limit of this hi- 
erarchy is obtained by replacing the differential opera- 
tor <p = hdt by a function ef> —* <f>, and the commutator 
ih~ 1 [, ] by the Poisson bracket { , }. We have seen that 
in this limit the formal objects of the Lax-Sato construc- 
tion acquire a clear physical interpretation in terms of 
the Hele-Shaw flow. Namely, cf>, L and M become respec- 
tively the complex potential </>, the coordinate z, and 
the height function h. This correspondence also gives a 
meaning of h as a quantum of area. 

10. Applications The integrable structure revealed in 
this letter allows one to use the KdV theory in the study 
of the Hele-Shaw and other moving boundary problems. 
Special solutions of the hierarchy describe (i) bubble coa- 
lescence, (ii) bubble break-off, (iii) branching, (iv) bubble 
creations, etc. We report some of them in |l4j . 

It is known that the dispersionless limit of non-linear 
waves is singular [l6| . The solutions suffer from nonphys- 
ical shock waves. Smooth initial data generally evolve 
into a multi- valued ( "overhanged" ) function within a fi- 
nite time. Such nonphysical solutions are equivalent to 
the finite-time singularities of the Hele-Shaw flow. 

The singular behaviour should be corrected by a regu- 
larization of the cusps at short distances, by introducing 
a new scale - a short distance cut-off. Surface tension, 
lattice or DLA-like regularizations were considered. 

The integrable structure of singularities suggests a 
novel, "dispersive" regularization. The flow is seen 
as an ill limit of the true dispersive flow, just as Eq. iQ 
is seen as a singular limit of Eq. I|14|l . We will present de- 
tails of this regularization elsewhere. Here we only note 
that this regularization treats the Hele-Shaw flow as a 
stochastic process of deposition of small particles with 
an area h. It carries a similar physics as the DLA |3j. 
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